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ABSTRACT 

We evaluate the dependence of the cluster correlation length r c on the 
mean intercluster separation D c , for three models with critical matter density, 
vanishing vacuum energy (A = 0) and COBE normalized: a tilted CDM 
(tCDM) model (n = 0.8) and two blue mixed models with two light massive 
neutrinos yielding tt h = 0.26 and 0.14 (MDM1 and MDM2, respectively). All 
models approach the observational value of <t 8 (and, henceforth, the observed 
cluster abundance) and are consistent with the observed abundance of Damped 
Lymana systems. Mixed models have a motivation in recent results of neutrino 
physics; they also agree with the observed value of the ratio og/ ^, yielding the 
spectral slope parameter T, and nicely fit LCRS reconstructed spectra. 

We use parallel AP3M simulations, performed in a wide box (side 
360 /i^Mpc) and with high mass and distance resolution, enabling us to build 
artificial samples of clusters, whose total number and mass range allow to cover 
the same D c interval inspected through APM and Abell cluster clustering data. 

We find that the tCDM model performs substantially better than n = 1 
critical density CDM models. Our main finding, however, is that mixed models 
provide a surprisingly good fit of cluster clustering data. 

PACS: 95.35; 98.80; 98.65. Cw 

Subject headings: dark matter: massive neutrinos, large scale structure of the 
universe, methods: numerical, galaxies: clustering, clusters 
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1. Introduction 

The study of the clustering of galaxy clusters, in the early eighties, allowed a basic 
advancement in our understanding of Large Scale Structure (LSS). The discrepancy between 
the galaxy correlation length r g and the cluster correlation length r c (Bahcall & Soneira 
1983, Klypin & Kopylov 1983, but see also Hauser & Peebles 1973) led to the introduction 
of the concept of bias (Kaiser 1984). Data on r c were then worked out, in further detail, for 
Abell clusters by Peacock & West (1992) and Postman, Huchra & Geller (1992), as well as 
for clusters in APM and in Edinburgh-Durham Southern Galaxy catalogs, by Dalton et al 
(1992), Nichol et al (1992) and Croft et al (1997). 

These analyses show that the value of r c depends on the mass threshold (M th ) of 
the cluster sample, through its mean intercluster separation D c = n _1 / 3 (> M th ), and 
that r c increases with D c . However, r c values obtained from Abell and APM data seem 
only partially consistent; this is to be partially ascribed to different cluster definitions; 
Bahcall & Burgett (1986), Bahcall & Cen (1992) and Bahcall & West (1992) suggested that 
observational ambiguities are wide enough to allow to conjecture that the scaling relation 
r c ~ 0.4 D c holds for 20 < D c h/Mpc < 100. Herebelow, we shall refer to this relation as 
BW conjecture. It ought to be born in mind that, above ~ 50/i _1 Mpc, such conjecture 
hinges on the estimates of r c for 55 and 94 /i _1 Mpc mean separations, for richness R > 1 
and R > 2 Abell clusters, while APM data, for the same D c range, give smaller r c . Dekel 
et al (1989) and Sutherland & Efstathiou (1991) suggested that the projection effects and 
peculiar inhomogeneities in the Abell sample might have biased upward r c at large D c . 
Peacock & West (1992), instead, confirmed such points (see also Jing, Plionis & Valdarnini, 
1992). Altogether, it may be fair to say that the controversy on the observational behaviour 
of r c for high D c values has not been solved yet, although, as we shall see, there may be 
good reasons to assess that different cluster definitions play a key role. 
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This paper is devoted to a comparison of cluster clustering, as it emerges from such 
data, with simulations of three cosmological models: a tilted CDM (tCDM) model and two 
mixed models (MDM1 and MDM2) with cold+hot DM. All models have critical matter 
density, vanishing vacuum energy, and are COBE normalized. During the last few years, 
much attention has been devoted to models with a positive cosmological constant A, also 
because of the remarkable data sets concerning SN la (see, e.g., Riess et al 1998, Perlmutter 
et al 1998, and references therein). In this work, we shall not debate whether mixed models 
can still offer a fair fit to all cosmological data; they certainly do not fit SN la data, unless 
their current interpretation was misled by some systematic bias. In a number of cases, 
however, mixed models were just not tested and the success of A-models to fit some data 
set was directly taken as further evidence in their favour. In the case of cluster clustering, 
we shall show that mixed models perform quite well and are surely better than any other 
model with matter density parameter Vt rn = 1 considered until now. 

In order to fit cluster data with a model, a large simulation volume is required; in fact, 
we need a fair sample of galaxy clusters for large D c , as well as adequate mass and force 
resolutions, to identify clusters in a reliable way, for small D c . Simulation parameters are 
therefore set so to allow a sample of 90 clusters, at least, for large D c and ~ 60 baryon-CDM 
particles per cluster, at least, for small D c (as we shall see, 60 particles correspond to 
~ 10 14 /i _1 M Q ). Altogether, at redshift z = 0, we shall therefore span a D c interval ranging 
from ~ 20 to 80 /i -1 Mpc. 

Cluster clustering has been studied by various authors in simulations. In particular, 
the behaviour of r c vs. D c , for standard CDM and open CDM, was studied by Bahcall & 
Cen (1992), Watanabe et al (1994), Croft & Efstathiou (1994), Eke et al (1996), Croft et 
al (1997), Governato et al (1999). Their results allow to conclude that CDM models with 
n — 1 may approach the observed behaviour of r c vs. D c , only for Vt m < 1. The behaviour 
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of r c vs. D c in a mixed model was also studied, using PM simulations, by Klypin & Rhee 

(1994) and Walter and Klypin (1996). Their work treated a different mix from those 
considered here, using smaller box and resolution. Accordingly, they could inspect only the 
D c interval running from ~ 20 to 45 h~ 1 Mpc. The behaviour they found is only marginally 
consistent with a constant r c /D c ratio, but their model does not exhibit much improvement 
in respect to pure CDM. 

The mixed models we consider here were selected on the basis of recent tests on v 
flavour mixing, which seem to support a non-vanishing z/-mass. Mixing data come from 
the solar v deficit (see, e.g., Hampel et al 1996, for GALLEX, and Abdurashitov et al 1996, 
for SAGE), the atmospheric v anomaly (Fukuda et al 1994) and the LSND experiment 
(Athanassopoulos et al 1995) on z/'s arising from fi + and ir + decay. Barger, Weiler & 
Whisnant (1998) and Sarkar (1999) show that all above results can agree if a fourth sterile 
v exists, which can be however added without harming BBNS or LEP standard results. 
Diagonalizing the mass matrix, they eventually obtain the four z/-mass eigenvalues, which 
split into two nearly degenerate pairs, corresponding to m u ~ and m v ~ 1.4-1.5 eV. It 
must be outlined that, within this picture, there remains no contradiction among different 
experimental results, at variance with earlier analyses which seemed to find contradictions 
between LSND and other //-mixing results. 

In a cosmo logical context, however, mixed models have been considered since long. The 
transfer function for several mixed models was first computed by Bonometto & Valdarnini 
(1984). Results on mixed models were then found by a number of authors (see, e.g., 
Achilli, Occhionero & Scaramella 1985, Valdarnini & Bonometto 1985, Holtzmann 1989, 
for results obtainable from the linear theory, and Davis et al 1992, Klypin et al 1993, 
Ghigna et al 1994, for early simulations). After the release of LSND data, Primack et al 

(1995) performed simulations of models with 2 massive z/'s and yielding fl h = 0.20 and 
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found that such mixture eased some problems met by greater Q h models. The possibility 
of considering mixed models together with blue spectra (primeval spectral index n > 1) 
was first considered by Liddle et al (1996) and Lucchin et al (1996). In the former 
paper, blue mixed models able to fit all linear and analytical constraints were shown to 
exist. In the latter paper, inflationary models leading to blue spectra were discussed and 
results of an N-body simulation of blue mixed models were reported. Unfortunately, the 
model considered violated some observational constraints. A systematic study of blue 
mixed models was recently performed by Bonometto & Pierpaoli (1998) and Pierpaoli & 
Bonometto (1999), selecting those consistent with CMB data and data predictable from the 
linear theory. 

In the next section we show that the models considered here, on the basis of //-physics 
motivation, are also suitable to fulfill the main observational constraints. In § 3 we review 
the technique used to simulate their non-linear evolution. In § 4 we describe how clusters 
are selected in simulations. Then, in § 5 we describe how the 2-point correlation function 
and its error estimates were worked out. In § 6, we will show the main results of the 
r c vs. D c behaviour derived from fits to the 2-point functions. § 7 is devoted to discussion 
of the results and the main conclusions we derived from this work. 



2. Model parameters 

The mixed models discussed in this paper were already considered in a previous work 
by Gardini, Bonometto & Murante (1999; hereafter Paper I) They are models with two 
equal-mass massive z/'s, selected by requiring agreement with data which can be fitted 
using the linear theory. More in detail, we required, first of all, agreement with observations 
at top and bottom scales, i.e. with COBE data and with the observed Damped Lymana 
system (DLAS) abundance (Storrie-Lombardi et al. 1995). Assuming 2 massive z/'s 
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with m u > 1.5 eV, we adjusted the spectral index n so to agree with the above top and 
bottom data, choosing the minimum allowed value for the spectral amplitude (Ay). Over 
intermediate scales, the main constraints to be tested are at 8 and 25 h~ 1 Mpc, where we 
evaluated the mass variances <T8,25 (see below). From such values we can work out the 
expected spectral slope and cluster abundance (again, see below); comparing their values 
with observations we see that m u values up to ~ 3 eV can be considered, without violating 
such constraints. Accordingly, we considered two values for the hot-dark-matter (HDM) 
density parameter Qh, yielding m v at the top and bottom of the allowed interval. 

A CDM model, selected so to fit the same data in a similar way, was also studied, for 
the sake of comparison. While mixed models require n > 1, even with low Ay, CDM may 
fit COBE data only if n< 1. 

Model parameters are shown in detail in Table I, while fig. 1 shows the spectra obtained 
from the linear transfer function T(k) against APM reconstructed spectral points ( |Baugh| 
fe Gatzahaga 1996]) . The wavenumber k is related to the comoving length scale L = 2ir/k 



and to the mass scale M = (47r/3)p L 3 , where p Q is the present density of the Universe. In 
Fig. 1 also other data and results are shown, which will be discussed below. 



EDITOR: PLACE TABLE HERE. 



The tCDM model approximates the APM galaxy spectrum slightly better than 
the standard CDM (sCDM), thanks to its increased slope. However, it lays still quite 
below the spectral points around the peak at k ~ 5 x 10~ 2 /iMpc _1 (comoving length 
A ~ 100-120 /i -1 Mpc). Blue mixed models, instead, show a stronger spectral peak, 
occurring where v free streaming bends a primeval spectrum steeper than Zel'dovich. This 
causes a steeper downward spectrum for 5 /i _1 Mpc^ 50 /i -1 Mpc, which might be related 
to the reasons why blue mixed models approach the r c vs. D c behaviour up to large D c . 
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EDITOR: PLACE FIGURE |T| HERE. 

Mass variances are defined according to the relation 

AL) = I ^A, fduu^T* (£) W\u) , (1.1) 
with a top-hat window function W(u) = 3(sinw — u cos u)/u 2 and 

2-7T 3 A t 

p{k) = —-£(*„*)" . (i.2) 

3 x A 

Here x Q is the co moving horizon distance. Using eq. (1.1) we estimate the parameter 

T = 7.13 x KTVsMs) 1073 , (1.3) 

which, only for pure CDM models, is often approximated as V ~ Q a h (see Efstathiou et 
al (1992). Peacock & Dodds (1994), using APM data, and Borgani et al. (1994) using 
Abell/ACO samples, constrained T within the (2a) intervals 0.19-0.27 and 0.18-0.25 
respectively. This parameter therefore tests the spectral slope above cluster mass scales. 

Values of cr 8 are directly related to the expected cluster number densities. A direct fit 
of data with simulations, for the cluster mass function, is given in Paper I, where, however, 
a different definition of cluster mass was used. This point will be discussed again below and 
will be deepened in a forthcoming work. 



3. The simulations 

The three simulations considered here were already used in Paper I, to which we refer 
for details. They were performed using a parallel N-body code, based upon the serial public 
AP3M code of Couchman (1991), extended in order to treat variable mass particle sets and 
used varying the time-steps, when needed. We considered a box of side L = 360 h~ 1 Mpc 
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(h is the Hubble parameters in units of 100 kms _1 Mpc _1 ); here CDM+baryons were 
represented by 180 3 particles, whose individual mass is m 180 = 2.22 x lO 12 /i _1 M for 
tCDM. Mixed models also involve 2 massive z/'s with m v ~ 3.02 eV and 1.63 eV, to yield 
Vl h = 0.26 and 0.14 (MDM1 and MDM2, respectively). Hence, slow particles, representing 
CDM+baryons, have masses f2 s m 180 with Q s = VL C dm + ^bar = 0.74 and 0.86, respectively, 
while fast particles, representing HDM, have masses ^Qh m iso (the ratio 2:1 between fast 
and slow particle numbers is required to set initial conditions with locally vanishing linear 
momentum). Our force resolution can be reported to a Plummer-equivalent smoothing 
parameter e p i ~ 40.6 /i _1 kpc. The comoving force and mass resolutions approach the limits 
of the computational resources of the machine we used (an HP Exemplar SPP2000 X 
Class processor of the CILEA consortium at Segrate-Milan) . The numerical resolution 
of our simulations were similar to other simulations of pure CDM, with different initial 
conditions performed by Colberg et al. (1997), Thomas et al. (1997), Cole et al. (1997), 
and Governato et al (1999). Mixed model simulations with a comparable dynamical range, 
instead, have only been performed by Gross et al. (1998), but in a smaller volume. 

4. Cluster selection 

Different criteria can be used to select clusters in simulations. In paper I, we identified 
clusters as virialized haloes. Here we shall give results obtained with a cluster definition 
aimed to approach more closely their observational definition. However, we widely tested 
and compared results ensuing different definitions, and a comprehensive discussion of the 
fit between particle sets obtained with various criteria will be published elsewhere. Let us 
however state that 2-point function outputs are fairly robust; although specific values of 
the clustering length r c , at various D c , have even significant variations, when the cluster 
definition is changed, the general trend is always preserved: the mixed models discussed 
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here however fit observational data. We shall illustrate this point with a few examples, 
without giving outputs for whole sets of different cluster definitions. 

In order to approach the observational pattern, here clusters were found with a 
spherical overdensity (SO) algorithm, based on a fixed sphere radius R a = 1.5/i _1 Mpc. 
The details of the procedure are close to those suggested by Croft & Efstathiou (1994) and 
Klypin & Rhee (1994). Hence, effects arising from the limiting magnitude of (observational) 
samples, border effects and projection effects are not included. (Also the sphere radius is 
fixed to mimic the Abell cluster definition; clusters found in the APM survey were selected 
also with smaller R a . Our simulations were also used to test whether systematic effects 
arise from a different choice of R a ; the differences we found have no significant impact on 
the results that will be shown here.) 

More in detail, we start the procedure with a FoF algorithm, which finds sets of N 
CDM-baryon particles closer than / times the average inter-particle separation. Results 
reported here are obtained using / = 0.2 and N = 25. Centers-of-mass (CM) of FoF groups 
are then inspected as possible centers for SO. Starting from them, we follow an iterative 
procedure: CDM-baryon particles within a distance R a from CM and their CM are found; 
this is repeated until we reach a stable particle set and fix their CM. Only particle sets 
containing at least 25 particles are however kept. When, during the iterative procedure, 
two spheres intersect, only the most massive particle set is kept. Our procedure aims to 
find all clusters above a suitable mass scale. Loose requirements were therefore set on /, 
in order to explore any possible matter condensation; the dependence of our results on N 
was also tested. Reducing N obviously leads to more FoF groups and a number of them 
survives the iteration procedure defined hereabove. Most of such extra-clusters, however, 
do not contain many particles. The result of such tests can be summarized by stating that 
extra-clusters of more than ~ 60 particles, found lowering IV down to 12, are less than 
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~ 0.3 %, in all cases; this percentage has no further increase when still lower N are taken. 
Henceforth, for N > 60, i.e., for M > 1.3f2 s lO 14 /i _1 M = M min , our cluster samples can be 
considered complete. In fig. 2, we show the relation between cluster masses and D c values, 
both at z = and at z — 0.8. 

EDITOR: PLACE FIGURE g HERE. 

Among other tests, we also verified the size of virialized haloes contained in clusters, 
as a function of their mass M c . Let R v be the radius encompassing a sphere, whereinside 
the density contrast is 180, found starting from the CM of each cluster, but whose actual 
center is attained through a suitable number of iterations, so to be the CM of all particles 
within R v from it. Let then M v be the mass of all CDM+baryon particles within R v . For 
large clusters, R v may exceed R a ; when M c is approximately < 7Q S ■ 10 1A h~ 1 M & , in general 
it is M v < M c . In fig. 3 we show the values of M v , as a function of M c , for tCDM. The 
trend is quite similar for mixed models. 

EDITOR: PLACE FIGURE | HERE. 

Although M v tends to increase with M c , the trend is clearly not monotonic; this is due 
to the spread of the values we find for R v at any M c . In spite of that, for M > M m i n , all 
clusters contain a virialized halo. However, if we order by mass the cluster set, using M v 
instead of M c , we find a different result. Hence, cluster sets, whose mean distance is D c , 
are different if we use M v instead of M c . It is then significant to compare the dependence 
of r c on D c for the two different orderings. In Paper I, clusters were given M v as mass; 
such definition is farther from observational criteria, but is likely to be closer to physical 
requirements, e.g., if we aim to compare simulation outputs with the expectations of a Press 
& Schechter approach. Herebelow, in a few cases, we shall test how the clustering length 
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depends on D c , when M v replaces M c . As we shall see, outputs depend significantly on the 
model, but are substantially independent from the cluster definition. 



5. Cluster 2— point correlation function 

Using clusters in our simulation box, ordered according to their M c , we computed the 
2-point correlation function £(r), for a set of D c values, by applying the estimator 

«M = - 1 • (3-D 

Here N pairs (r) is the number of cluster pairs in the radial bin of volume SV(r), centered 
on r, and L 3 is the box volume. Error bars for £(r) were estimated using the standard 
bootstrapping procedure. (We checked the convergence of the estimator of the standard 
deviation evaluated from bootstrap realizations, by inspecting the third moment of the 
bootstrap distribution. In all cases, convergence was attained when the number of bootstrap 
realizations matched the number of points in the catalogues; see e.g. [Bradley 1982j ) We 



compared such errors with the usual Poisson errors, which were found to be systematically 
smaller by a factor ~ 2. 

We then performed two different fits to a power law 

£(r) = (r/r c n , (3.2) 

over the distance range 5/i _1 Mpc< r < 25/i -1 Mpc: (i) A constrained fit, assuming a 
constant 7 = 1.8. For the sake of example, in fig. 4 we show such fit for D c = 30 /i -1 Mpc. 

EDITOR: PLACE FIGURE fH HERE. 



(ii) An unconstrained fit, allowing both r c and 7 to vary. Points were weighted by 
the corresponding bootstrap errors and r c best-fit values are also given with bootstrap 
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errors. Such errors are obviously smaller for the constrained fit, where our ignorance on 7 is 
hidden. In general, large D c clusters yield best fit 7 values approaching 2, although 1.8 lays 
always within la. Our r c estimates are performed at z = and z = 0.8, to inspect cluster 
clustering evolution. 

6. Results 

In fig. 5 we report the r c vs. D c behaviour for tCDM, MDM1 and MDM2, for fixed 
7 = 1.8. In fig. 6 we give results for the same cases, obtained with 2-parameters fits on r c 
and 7. Errors bars represent 1 a bootstrap errors (see § 3). Of course, error bars are smaller 
in the single parameter fits, where our ignorance on 7 is hidden. 

Together with the r c values obtained from our simulations we also plot APM and Abell 
cluster data, the BW conjecture, and the results from simulations performed by Bahcall 
& Cen (1992) and Croft & Efstathiou (1994). Recent results obtained by Governato et al 
(1999), for a critical CDM model, lay between the last two curves. Observational points 
and error bars given in our figures were obtained from original work. We draw the reader's 
attention on the fact that, in some recent work studying cluster clustering in simulations, 
observational points and error bars are not accurately reported. 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE | HERE. 

A comparison of tCDM with data, shows that simulated and APM data points are 
in fair agreement. In view of the better fits obtainable with mixed models, shown in the 
same figures, one might tend to overlook the improvement of tCDM in respect to CDM 
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models with n — 1, which, instead, is significant. Our tCDM model, however, seems to miss 
systematically the Abell catalog points and, thence, is far from the BW conjecture, which 
tries to set a compromise between Abell and APM results. 

From this point of view, the performance of mixed models is better. For low D c , 
MDM1 tends to give r c values above the BW conjecture (see, however, Lee & Park 1999). 
A similar, but less pronounced, effect exists also for MDM2. On intermediate scales MDM1 
sticks on the BW conjecture curve and meets two of the APM points at the 2-cr level only. 
MDM2, instead, seems to try to compromise between APM and Abell points. On top 
scales, the behaviours shown here by the two models are opposite. The MDM1 behaviour, 
at such scales, seems however somehow anomalous; such scales are those which are most 
likely affected by cosmic variance and the MDM1 behaviour at D c > 65-70 ft, _1 Mpc should 
certainly be tested with different model realizations. Furthermore, unconstrained fits tend 
to indicate that such discrepancy arises from different correlation function slopes. 

It may also be significant to consider the unconstrained fit obtained ordering clusters 
according to M v masses, which is shown in fig. 7. Let us notice that: (i) in most cases, 
error bars are smaller; (ii) the peculiar feature for MDM1 at large D c has disappeared. It is 
likely that such improvement is related to a more direct physical significance of the mass 
M v and the (variable) radius R v . Taking a fixed radius R a , instead, risks to accentuate a 
dependence on local peculiar features. In principle, this is more likely to occur for small 
mass clusters, for which significant volumes, still unaffected by virialization processes, lay 
within R a . In our simulation volume, however, we have a large deal of low-mass clusters 
and this allows an efficient averaging over local realizations. At the top mass end, instead, 
the sample is more restricted and we must mostly rely on the virialization process, rather 
than on sample averaging, to smear off local peculiarities. Our results seem to indicate that 
significant memory of initial conditions is kept also below R v . 
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EDITOR: PLACE FIGURE Q HERE. 

In Fig.s [8] and |9| we report a comparison between the 2-point function results at z = 
and z = 0.8, obtained using M c . 

EDITOR: PLACE FIGURE |8| HERE. 



EDITOR: PLACE FIGURE | HERE. 

Unconstrained fits at z = 0.8 are rather noisy at large D c , in particular for top scales. 
Constrained fits, instead, might be taken as an indication of clustering evolution on top 
scales. Here, perhaps, there is a further evidence of anomaly in MDM1, which is the only 
case when clustering seems however weeker at z = than at z — 0.8. Apparently, all models 
seem to indicate a greater clustering length at scales between 50 and 65-70 /i _1 Mpc for 
z = 0.8. Apart of MDM1, instead, this is inverted above 70 /i -1 Mpc. 



EDITOR: PLACE FIGURE |T0| HERE. 



Fig [L(], instead, shows results on cluster evolution based on M v ordering and using 
constrained fits. The kind of evolution found hereabove seems confirmed, while MDM1 
anomaly is reduced. 

7. Discussion 

Previous numerical results on cluster clustering, based on models with Q m = 1, gave a 
behaviour of r c vs. D c a few it's below observational results. The only exception are Bahcall 
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& Cen (1992), whose numerical study involves peculiar extrapolations and however succeeds 
to meet two APM points at large D c only. When considering subcritical CDM models 
(OCDM), the same authors obtain a behaviour close to the BW conjecture. However, 
this is not fully confirmed by later numerical studies; although they clearly indicate that, 
in OCDM models, the D c dependence on r c is consistent with APM points, Abell cluster 
points seem to require a still steeper dependence than in OCDM. Such findings, however, 
were currently interpreted as an indication that observational data on cluster clustering 
may be approached only by models with Q m < 1 and led to arguing that the observed 
dependence of r c on M t h is somehow related to an early cluster formation. 

The behaviour we find for tCDM does not support such kind of inference. Taking a 
spectral index n ^ 1 has no substantial effect on the time of cluster formation, which is 
quite similar to standard CDM. In view of the more stricking outputs for mixed models, we 
must not disregard the result we find for tCDM. The r c vs. D c behaviour of clusters in such 
model is analogous to previous outputs for OCDM and lays well above the very behaviour 
obtained by Bahcall & Cen (1992) for standard CDM. Our findings are that, taking n < 1, 
the clustering of clusters in an Vt m = 1 model approaches the behaviour obtained from APM 
cluster data. 

Even more striking is the cluster clustering behaviour for the mixed models considered 
in this work, which are based on low mass z/'s. In this case, the slope of the r c vs. D c 
behaviour approaches the BW conjecture and, more significantly, within 1 or 2-cr bootstrap 
error bars, we mostly find consistency both with APM and Abell results. The only exception 
could be the point of Abell clusters of richness R > 2, which is approached only by MDM2. 
Our results, however, support previous claims that, at the top mass end, wider samples 
may be required to suppress the cosmic variance. 

Quite in general, it can be reasonable to consider cluster clustering as a measure of 
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the spectral power on scales exceeding ~ 25 h^Mpc. All models with Vl, m = 1 ought to 
have similar values of a 8 , in order to be consistent with the observed cluster abundance. 
Accordingly, the power at scales ~ 25 /i~ 1 Mpc is basically gauged by the value of the T 
parameter. It may not be a case that models perform in a better way as their T's approach 
the observational interval. 

Surely, in the case of mixed models, the situation is complicated by the presence of a 
hot component, which may slow down the gravitational growth in the non-linear regime, 
in a scale-dependent fashion. Previous results for mixed models concerned a mix including 
30 % of HDM, due to a single v with mass m u ~ 7eV. Here we deal with z/'s 3-4 times 
lighter; accordingly, when cluster formation begins, their speeds are ~ 3-4 times greater. 
Peculiar effects of MDM are therefore significantly reinforced. Hence, besides of having a 
suitable spectral slope, the mixed models treated here are still more different from standard 
CDM, because of the late v derelativisation. 

In this paper we showed that critical CDM models, with a blue spectrum suitably 
"compensated" by a light— i/ component, besides of fitting most LSS and CMB data, are 
able to follow APM and Abell cluster clustering data. This result adds to those discussed 
in Paper I, where critical blue mixed models were shown to provide a good fit to the cluster 
mass function and to be in agreement with Donahue at al. (1998) findings concerning 
high-z cluster abundance. 
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Table 1: Parameters of the models; besides of input parameters or quantities derived 
from the linear theory we report the value of N c i (sim) (number of clusters of mass 
> 4.2 x 1O 14 /i~ 1 M in a box side of 100/i _1 Mpc) as obtained in Paper I. The normalization 
to COBE quadrupole was deliberately kept at the ~ 3 a lower limit, in order leave some 
room to the contribution of tensor modes, but keeping however consistent with data. The 
expected interval for N d is 4-6, but models with N d up to 8-10 cannot be safely rejected. 
L a = £l b £l co n(z = 4.25, M = 5 x 10 9 /i _1 M Q ) x 10 3 accounts for the amount of gas expected in 
Damped Lyman-a systems. More details can be found in Paper I. Data provided by Storrie- 
Lombardi et al. (1995) give L a > 2.2 ± 0.6. This is one of the most stringent indicators of 
the model capacity to produce high-z objects. 
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Fig. 1. — Spectra of the three models at z — 0. Solid curves give the linear power spectrum 
and the spectrum corrected for non-linearity, according to Peacock & Dodds (1996). Empty 
squares yield the simulation spectra corrected for CIC (cloud-in-cell; see paper I for more 
detail). Circles with 2 a errorbars are the power spectrum measured from the APM survey. 
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Fig. 2. — Mass limit (M th ) of cluster samples with mean separation D c = n 1 / 3 (> M ; 
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Fig. 3. — All clusters with M c > 10 l4 h 1 M© are found to contain a virialized halo. Its mass 
My is shown, as a function of M c , for the tCDM model. 
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Fig. 4. — The 2-point correlation function estimate (solid points) and la bootstrap error 
bars for D c = 30 h~ x Mpc. A constrained fit with 7 = 1.8 is also shown (solid line). 
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Fig. 5. — The cluster correlation length r c as a function of D c , for the three models. Here 
r c is obtained by requiring 7 = 1.8. la bootstrap error bars are also shown. r c values 
obtained from APM and Abell cluster data are also reported (filled circles and triangles, 
respectively). The thick solid line corresponds to the BW conjecture. Results of SCDM 
simulation of Bahcall & Cen (1992) are shown by the thin solid line, while dashed and 
dotted lines refer to Croft et al (1997) simulations of LCMD and SCDM models. 
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Fig. 6. — Same as Fig |5] but with r c obtained by simultaneously fitting it and 7 to the halo 
distribution. Error bars, representing la deviation derived from the weighted least-square 
fits, are obviously larger than in the fixed-7 case. 
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Fig. 7. — Same as Fig |] but with D c obtained ordering clusters according to M, 
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Fig. 8. — Comoving clustering evolution from z = 0.8 to z = obtained by fitting r c with 
7 = 1.8. 
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Fig. 9. — Same as Fig|8|, but from simultaneous fits of r c and 7. MDM1 shows some anomaly. 
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Fig. 10. — The same as Fig|8|, but D c are obtained ordering clusters according to M v . Plots 
are less noisy, when such mass determination is used. 



